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^h ■ Abstract 

An increasing number of studies use the spectrum of cardiac signals for an- 
alyzing the spatiotemporal dynamics of complex cardiac arrhythmias. How- 
ever, the relationship between the spectrum of cardiac signals and the spa- 
tiotemporal dynamics of the underlying cardiac sources remains to date un- 
clear. In this paper, we derive a mathematical expression relating the spec- 
trum of cardiac signals to the spatiotemporal dynamics of cardiac sources 
and the measurement characteristics of the lead systems. Then, by using 

t-h \ analytical methods and computer simulations we analyze the spectrum of 

cardiac signals measured by idealized lead systems during correlated and 

\q ■ uncorrelated spatiotemporal dynamics. Our results show that lead systems 

can have distorting effects on the spectral envelope of cardiac signals, which 
depend on the spatial resolution of the lead systems and on the degree of 
spatiotemporal correlation of the underlying cardiac sources. In addition to 
this, our results indicate that the spectral features that do not depend on the 
spectral envelope, such as the dominant frequency, behave robustly against 
different choices of lead systems. 
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1. Introduction 

Fibrillation is a complex cardiac arrhythmia whose spatiotemporal char- 
acteristics remain poorly understood. Traditionally fibrillation has been de- 
scribed as random and disorganized, since it induces highly irregular traces in 
the electrocardiogram (ECG) signal. However, with the application of non- 
linear dynamics theory to the investigation of cardiac arrhythmias, and the 
development of optical and electrical mapping techniques, it has been sug- 
ested that fibrillation can possess some degree of spatiotemporal regularity 
l|, [2] . This view has led in a natural way to study fibrillation based on the 
spectrum of cardiac signals such as the ECG and intracardiac electrograms 
(EGM). Spectral features of cardiac signals have been proposed as experi- 
mental indices for detecting ventricular fibrillation (VF) [3j, for quantifying 
the degree of spatiotemporal organization of atrial fibrillation (AF) |4| and 
for predicting the success of defibrillation shocks p, la, LZ|. Also, intracardiac 
mapping techniques have been combined with dominant frequency (DF) anal- 
ysis to study the spatiotemporal characteristics of fibrillation. This method, 
known as DF mapping, has revealed spatiotemporal regularities during AF 
in both animal models |S|, |9|, [ly] and in patients 111, ll2[ and it is currently 



regarded as a potential technique to guide AF ablation therapies [13 . 

Despite the increasing number of studies that use spectral techniques to 
analyze fibrillation, the meaning of the spectrum of cardiac signals remains 
to date elusive. Even though individual spectral features of cardiac signals 



have been linked to spatiotemporal characteristics of cardiac rhythms [14 . 



15l . [la], the relationship between the spectrum of cardiac signals and the 
spatiotemporal characteristics of cardiac rhythms has not been thoroughly 
investigated. In addition to this, the effects of lead systems on the spectrum 
of cardiac signals are not well understood, and consequently it is not clear 
how the spectra of cardiac signals measured by different lead systems relate 
to one another. The elucidation of the relationship between the spectra 
of cardiac signals measured by different lead systems is of technical and 
clinical interest in the context of fibrillation, since it would contribute to the 
development of novel, improved methods of DF cardiac mapping, such as 



non-contact intracardiac electrical mapping [121, [18] and non-invasive surface 
ECG mapping 0. 

In this paper, we develop a mathematical formalism for investigating, in 



a systematic way, the spectral manifestation of different cardiac rhythms and 
the spectral effects of the lead systems. By following a multivariate signal 
analysis approach, we identify the mathematical relationship between the 
spectrum of cardiac signals and the spatiotemporal dynamics of the under- 
lying cardiac rhythms. Our formalism allows us to derive theoretical results 
which are relevant for the analysis and interpretation of the spectrum of 
cardiac signals, and for devising spectral methods for investigating the spa- 
tiotemporal dynamics of cardiac rhythms. 

The organization of this paper is as follows. In Section [2] we develop our 
mathematical formalism and connect the spectrum of cardiac signals to the 
spatiotemporal dynamics of the underlying cardiac rhythms. Physiologically 
meaningful cases are studied analytically in Section [3J and in a computer sim- 
ulation environment in Section HJ Finally, Section [5] conveys the conclusions 
of our investigation and the discussion. 

2. Mathematical formalism 

In this section we present the mathematical formalism for investigating 
the spectrum of cardiac signals. Firstly, we introduce the lead-field bioelec- 
tric model of cardiac sources and signals. Then, based on multivariate signal 
analysis, we define the autocorrelation and the spectrum of cardiac sources. 
Finally, we identify the relationship between: the spectrum of cardiac sig- 
nals, the spatiotemporal dynamics of cardiac sources and the measurement 
characteristics of lead systems. 

Throughout this paper the following notation is used: (-)t denotes time- 
average, J-[-] is the Fourier Transform operator, (*) denotes convolution and 
5(-) is the Dirac's delta. We use the following vector definitions: 1 = [1, 1, 1] T 
andO= [0,0, 0] T . 

2.1. Bioelectric model 

Cardiac sources are the bioelectric processes generated by the heart dur- 
ing contraction. There exist different, equivalent mathematical paradigms 
to model the activity of cardiac sources, such as the monopole field and 



the dipole field [20j. In this study, we model cardiac sources as a time- 
varying dipole field, i.e. as a spatial distribution of time-varying dipoles 
3(v,t) = [J x (v,t), J y (v,t), J z (v,t)] T on a volume V, where v denotes a point 
located inside V and t denotes the time instant. 



The time-varying activity of cardiac sources can be measured by lead 
systems, producing cardiac signals. Taking the dipole field as our reference 
description for cardiac sources, we follow a lead-field approach to model car- 
diac signals |20|. According to the lead- field theory, the cardiac signal c(t) 
that is induced at a lead system by a dipole field J(v, t) can be expressed as 

c{t) = f L T (v)3(v,t)dv (1) 

Jv 

where the vector field L(i>) = [L x (v),L y (v),L z (v)] T is the measurement sen- 
sitivity distribution (MSD) and describes the ability of the lead system to 
measure cardiac dipoles located at v G V. In words, cardiac signals are a 
weighted linear combination of the underlying cardiac sources. 

2.2. Autocorrelation and spectrum of cardiac sources 

The autocorrelation of a cardiac source, p(v, w, r), Wv , w G V, is defined 
as the collection of the cross-correlations between all pairs of dipoles in V. 
Since cardiac dipoles are vectorial entities, the cross-correlation between two 
dipoles consists of the cross-correlations between all three components of each 



dipole [2lj. In order to define the autocorrelation of a cardiac source, the 



average dipole field J(i>) needs to be introduced: 

3{v) = (J(v,t)) t 

= [{J x (v,t)) t ,(J y (v,t)) t ,{J z (v,t)) t ] T . (2) 

Based on J(v), we define the zero-average dipole field J'(v, t) = J(v, t) — J(t> ), 
so that (J'(v,t)) t = 0. The cross-correlation matrix between two cardiac 
dipoles J(f , t) and J(u>, t), where v , w G V, is then defined as 

p(v,w,r) = (J'(v,t + T)3' T (w,t)) t 

p xx (v,w,r) p xy (v,w,r) p xz (v,w,r) \ 

Pyx(v,w,T) p yy (v,w,r) p yz (v,w,r) . (3) 

P zx (v,w,t) p zy (v,w,r) p zz (v,w,T) J 

Therefore, each entry of p(v, w, r) contains the cross-correlation between one 
component of 3(v , t) and one component of J(u>, t). For instance, matrix en- 
try p zy (v,w,r) is (J' z (v,t + r)Jy(w,t))t- Also, the average power of dipole 
component J' x (v , t) is by definition P x (v) = p X x(v, v, 0), and analogous expres- 
sions can be obtained for the average power of dipole components J y (v, t) and 
J' z (v,t). 



The spectrum of a cardiac source, cr(v,w,f), Wv,w G V, corresponds to 
the collection of the cross-spectra between all pairs of dipoles in V, and is 
defined as 

<r(v,w,f) = F\p{v,w,r)\ 

a xx (v,w,f) a xy (v,w,f) a xz (v,w,f) \ 

a yx (v,w,f) cr yy (v,w,f) <j yz (v,w,f) (4) 

a zx (v,wj) a zy (v,w,f) a zz (v,wj) J 

where the operator J-[-] is applied to p(v, w, r) on a component-by-component 
basis. For instance, a zy (v,w, f) is J-[p zy (v, w,r)]. 

We also define the total cross- correlation Rj(v,w,r) between two car- 
diac dipoles 3(v,t) and 3(w,t) as the sum of the entries of p(v,w,r) and 
the total cross- spectrum Sj(v,w, f) as the Fourier Transform of Rj(v, w, r). 
Mathematically, they can be expressed as 

Rj(v,w,t) = l T p(v,w,r)l, (5) 

Sj(v,wJ) = l T a(v,wJ)l. (6) 

Finally, we define the normalized cross-correlation p(v, w, r) as the matrix 
of entries 

~ / x p ab {v,W,T) 

Pab (v, w, t) = , (7) 

where a,b E {x, y, z}, and the normalized total cross-correlation Rj(v, w, r) 
as 

Rj{v, W, t) = —— rr. (8) 

ma.x{Rj(v } w, t)\ 

T 

2. 3. Autocorrelation and spectrum of cardiac signals 

Let c(t) be a cardiac signal measured by applying L(t>) to a cardiac source 
of autocorrelation p(v,w,r) and spectrum a(v,w,f) in V. Define c'{t) as 
the cardiac signal c(t) minus its time-average value c = (c(t))t, 

c\t) = c(t) - c. (9) 

The autocorrelation function R c (t) of the cardiac signal c(t) is defined as the 
following average 21] : 



R c (r) = (c'(t + T)c / (t)) t , (10) 



and its power spectrum S c (f) is denned as the Fourier Transform of its 
autocorrelation function, 

S c (f) = F[R c (t)}. (11) 

Based on flTJ, the following relationships can be derived between R c (t) and 
p(v,w,r), and betwen S c (f) and cr(v,w,f) (see Appendix A): 

R c (t) = / L T (v)p(v,w,r)L(w)dvdw, (12) 

JVxV 

S c {f) = f \7{v)(r{v,wJ)~L{w)dvdw. (13) 

JVxV 

Equations (TT21) and (IT3"|) reflect the linear relationship between cardiac sig- 
nals and sources [c.f. flTJ], and can be used to gain insight into the nature 
of the autocorrelation and spectrum of cardiac signals. According to (J13)) . 
two factors determine the spectrum of cardiac signals. The first factor is the 
spectrum cr(v,w,f) of cardiac sources. It is worth noting that this is the 
solely feature of the spatiotemporal dynamics of cardiac sources that mani- 
fests on the spectrum of cardiac signals. The second factor is the MSD of the 
lead system, L(t>). Since L(i>) is specific for each lead system, (fT3l) reveals 
that cardiac signals measured by different lead systems will in general have 
different spectra for the same underlying spatiotemporal dynamics. 

3. Analytical study 

In this section we study analytically the spectral manifestation of two 
second-order models of cardiac sources, namely the fully correlated (FC) 
source and the fully uncorrelated (FU) source. The FC and FU models 
are physiologically meaningful and can be used to describe the dynamics of, 
respectively, highly organized and highly disorganized cardiac rhythms. 

We firstly define the autocorrelation and the spectrum of the following 
models of spatiotemporal dynamics: identically distributed (ID), FC and FU. 
The ID model is introduced for facilitating the comparison of the spectrum of 
cardiac signals measured during FC and FU dynamics. Secondly, we define 
a simple, idealized model of MSD, namely the pulse model. Because of its 
simplicity, the pulse model is used in the analytical derivations and in the 
simulation experiments throughout this study. Finally, we derive analytically 
the spectrum of cardiac signals measured by pulse MSD during FU and FC 
dynamics. 
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3.1. Models of spatiotemporal dynamics 

We subsequently present the three second-order models of spatiotemporal 
dynamics that we use in this study, namely the ID, the FC and the FU 
models. 

3.1.1. Identically distributed spatiotemporal dynamics 

In this model of spatiotemporal dynamics all cardiac dipoles have the 
same autocorrelation and spectrum, 

p (v,v,t) = p{r), (14) 

<r(v,v,f) = *(f). (15) 

By substituting (JTlj) and (IT51) respectively into (J5J) and (jSJ), it can be proved 
that the total autocorrelation and total spectrum of all the dipoles are also 
identical, 

Rj(v,v,t) = Rj(r) = l T p(r)l, (16) 

Sj(v,vJ) = Sj{f) = l T a{f)l. (17) 

Note that this model only describes the activity of cardiac dipoles individu- 
ally, and does not specify p(v, w, r) nor <r(v, w, f) for v ^ w. 

3.1.2. Fully correlated spatiotemporal dynamics 

This model of spatiotemporal dynamics corresponds to highly regular 
rhythms, such as sinus rhythm, in which the activity of one dipole 3(w,t) 
can be expressed as a delayed version of the activity of another dipole 3(v,t), 

J(w,t) = J{v,t-{(v,w)), (18) 

where ((v, w) is defined as the time delay between the activities of dipoles 
3(v,t) and J(w,t). Based on (TT8|) . it can be proved (see Appendix B) that 
FC dynamics are also ID, p(v,v,r) = p{j) and cr(v,v,f) = cr(f) [cf. (fill) 
and (TT5"j) ]. and that the autocorrelation and the spectrum of FC sources can 
be expressed as: 

p(v,w,r) = p(r-((v,w)), (19) 

<r{v,w,f) = <r(f)exp[-j2irfC(v,w)]. (20) 

Consequently, FC sources are completely characterized by p(r), cr(f) and 



3.1.3. Fully uncorrelated spatiotemporal dynamics 

This model of spatiotemporal dynamics constitutes an idealization of 
highly irregular and disorganized rhythms, such as fibrillation, in which there 
is no second-order relationship between the temporal activity of any pair of 
cardiac dipoles. The autocorrelation and the spectrum of a FU source are 
defined as 

p(v,w,r) — p(v,v,r)S(v — w), (21) 

<t(v,w,t) — ct(v,v,t)S(v — w). (22) 

In words, the cross-correlation between two cardiac dipoles 3(w, t) and 3(v,t), 
where v ^ w, is null. 

3.2. Pulse model of MSD 

In this section we define one simple, idealized model of MSD, namely the 
pulse model. The pulse model describes a lead system that measures with 
the same sensitivity every dipole within a region Vq of the volume source V, 
while rejecting the rest. Mathematically, this model is defined as 

t f \ / l if u e VJ, 

L v {v) = < . ■ (23) 

I otherwise 

The pulse model can be treated as an approximation of physical MSD that 

effectively concentrate their measurement in a region Vq. 

For the subsequent analysis it is also convenient to quantify the spatial 

resolution (SR) of pulse leads. The SR can be defined as the region of the 

cardiac source that contributes the most to the measured signal. In this 

study, we quantify the SR of pulse leads by introducing the notion of the 

lead equivalent volume (LEV), which is defined as the relative size of Vo to 

the size of V, 

L dv Mm 
LEV = ^— = ^, (24) 

J v dv M v ' K J 

where My and My are the sizes of Vo and V respectively. Thus, for local 
measurements the LEV is close to zero, whereas for global measurements 
where Vq — V, the LEV is close to one. 



3.3. Spectrum of cardiac signals during FC and FU dynamics 
3.3.1. Fully correlated spatiotemporal dynamics 

We subsequently derive the spectrum of cardiac signals measured by pulse 
leads during FC dynamics. By substituting ( 123]) and (1201) into (1T3]) and using 
( 117]) . we can express the spectrum S c (f) in terms of the total spectrum Sj(f) 
and the time-delay ((v,w): 

S c (f) = / L^ (v)cr(v,w,f)L Vo (w)dvdw 

JVxV 

l T cr(f) exp[— j2Trf((v, w)]ldvdw 

■Vo 

= l T cr(f)l / exp[—j2irf((v,w)]dvdw 

JVqxVo 

= Sj(f) f exp[-j27rfC(v,w)]dvdw. (25) 

JV xV 

We can integrate (125]) with respect to ( by introducing the time-delay density 
function (TDDF) over Vo, F Vo ((). The TDDF F Vo (Q describes the frequency 
that a time-delay C is observed between two dipoles in Vo, and can be cal- 
culated as follows. Consider all the pairs of dipoles that reside in Vo- Given 
a time delay d, select all the pairs of dipoles whose activities are delayed 
by Ci- Then, the quantity -FV (Ci) corresponds to the density (fraction) of 
those pairs of dipoles, with respect to all the pairs of dipoles in Vo- Based on 
Fy (C), we can express (1251) as 

/oo 
exp[-j2nfC}F Vo (C)dC 
-oo 

= Sj(f)F[F Vi) (0} 

= Sj(f)S Vo (f). (26) 

In general, time delays ( are expected to be short for small Vo, i.e. small 
LEV [c.f. (1211) ]. Hence, the smaller Vo, the more concentrated -^(C) around 
( = and the broader Sv (f)- Therefore (l2"fiT) reveals that cardiac signals 
from local measurements (small LEV) during FC dynamics have broad band- 
widths, whereas cardiac signals from global measurements (large LEV) have 
narrow bandwidths. Equation (1261) is a natural result, since the integral in 
([1]) averages delayed signals at different points within the effective measure- 
ment region Vo- As a result, when the LEV is large the lead system acts as 
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a low-pass filter. By contrast, when Vq is small little averaging is performed 
and the distortion introduced by the lead system is low. 

3.3.2. Fully uncorrelated spatiotemporal dynamics 

We now derive the spectrum of cardiac signals measured by pulse leads 
during FU dynamics. By substituting ( 122|) and (123|) into ( 1T3|) . and using (|6j), 
the spectrum can be written as 

S c (f) = / L^ Q (v)cr(v,w,f)L Vo (w)dvdw 

JVxV 

= [ l T a(v,v,f)ldv 
JVo 

Sj(v,v,f)dv. (27) 

V 

In words, the spectrum of cardiac signals measured during FU dynamics is 
an average of the total spectra of the dipoles contained within Vq. If the FU 
dynamics is also assumed to be ID (ID-FU), then Sj(v,v,f) = Sj(f) [cf. 
([T7D1 and (ETD reduces to 



S c (f) = / Sj(f)dv = Sj(f)M Vo (28) 

JV 

Hence, except for the factor of My , the spectra of cardiac signals measured 
by pulse leads during ID-FU dynamics are identical to the total spectrum 
Sj(f). 

4. Simulation experiments 

In this section we present the results obtained in a computer simulation 
environment, designed for investigating the spectral manifestation of cardiac 
rhythms and the spectral effects of lead systems. Firstly, we describe the fol- 
lowing aspects of our simulation experiments: the cardiac source model, how 
we simulated FC and FU dynamics, and how we synthesized cardiac signals. 
Secondly, we use multivariate signal analysis to describe the simulated FC 
and FU dynamics. Thirdly, we analyze the TDDF of the simulated FC dy- 
namics and its relationship with the spectrum of measured cardiac signals. 
Finally, we analyze the effects of different pulse leads on the spectrum of 
cardiac signals during FC and FU simulated dynamics. 
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4-1. Simulation set-up 

Cardiac activity was modeled by following a probabilistic cellular au- 
tomata approach proposed in [22J which allows simulating complex spa- 
tiotemporal dynamics such as fibrillatory conduction and rotors. In this 
model, cardiac tissue is described as a grid of connected, excitable cells, 
characterized by a time-varying transmembrane voltage. Cell excitation pro- 
duces a transient change in the transmembrane voltage, causing excitation 
in the neighboring cells. This activity results in excitation waves traveling 
through the cardiac tissue. 

Based on this model, a 2D square sample of cardiac tissue consisting of 
101 x 101 cells was defined and two types of spatiotemporal dynamics were 
simulated. The first type of dynamics was generated by stimulating one side 
of the tissue sample at a rate of 1 Hz. This resulted in a FC dynamics 
consisting of a regular train of plane excitation waves traveling along the x 
axis (Fig. [1] (a)). The second type of dynamics corresponded to a chaotic 
excitation pattern characteristic of FU sources. This spatiotemporal dynam- 
ics consisted of multiple, fragmented excitation waves, traveling randomly in 
every direction (Fig. [1(b)). 

[Figure 1 about here.] 

By computing the voltage difference between neighboring cells, a 2D time- 
varying distribution of cardiac dipoles J(v,t) = [J x (v,t), J y (v,t)] T was sim- 
ulated. The amplitude of the simulated cardiac dipoles was highest at the 
wavefronts, where the voltage gradient is largest, and their direction was as 
expected aligned with the direction of propagation of excitation waves. In 
order to synthesize cardiac signals, eight square pulse MSD A n , 1 < n < 8, 
of increasing size and centered at the middle of the tissue sample were imple- 
mented. Table [T] shows the sizes of each pulse MSD along with their LEV, 
and Fig. [2] (b) depicts for illustrative purposes pulses A 5 and A 8 . Finally, 
a set of cardiac signals c n {t) were synthesized by applying each pulse A n to 
the simulated time-varying dipole distributions. 

[Table 1 about here.] 

4-2. Analysis of the spatiotemporal dynamics 

In order to analyze the spatiotemporal characteristics of both simulated 
dynamics, we obtained the average power of each dipole component, together 
with the normalized cross-correlation [c.f (J7|)] and the normalized total cross- 
correlation [c.f (JHJ)] between every pair of dipoles. 

11 



4-2.1. Fully correlated spatiotemporal dynamics 

In this dynamics, the normalized autocorrelation and total autocorre- 
lation were found to be approximately the same for every cardiac dipole, 
p(v,v,r) ~ p(r) and Rj(v,v,t) ~ Rj{t). Accordingly, the dynamics were 
classified as ID [c.f. ( TT4"|) and f lTB]) ]. Also, the average power P x {v) was found 
to be twice as high as P y {v). This observation is in agreement with the fact 
that most of the dipoles were aligned with the direction of propagation of 
the excitation waves, namely the x axis. 

Figure [1] (c) shows the normalized cross-correlations between the x and 
y components of dipole Vi, and dipoles V\ (itself), v 2 and v 3 . As expected, 
Pxx{ v i> v ii T ) revealed a high degree of periodicity, which could be approxi- 
mated by a regular train of impulses of period 1 s. This observation agrees 
with the simulated stimulation rate of 1 Hz. By contrast, p yy (vi,vx,T) 
showed a lower degree of periodicity since J y (v,t) components were ac- 
tivated in a more sporadical manner. The normalized cross-correlations 
Pxx(vi,v 2 , t) and p xx (vi,v 3 ,r) also consisted of a regular train of impulses 
of period 1 s and additionally, their maximum values occured at r = 0.004 
s and r = 0.048 s respectively. These delays corresponded to the travel- 
ing times of the excitation waves from v\ to v 2 and from v\ to v%, respec- 
tively. Thus, Pxx{vi,V2,t) ~ PxxiviiV^T - 0.004) ~ p xx {r - 0.004) and 
Pxx(vi,v 3 ,r) ~ p xx (vx, v 1 ,r - 0.048) ~ p xx (r - 0.048). In general, the maxi- 
mum value of the cross-correlation between any two dipoles occurred at a time 
delay ( that corresponded to the traveling time of the excitation wave from 
one dipole to the other, and also p xx (v, w, r) ~ p xx (v, v,t — £) = p xx {j — (). 
This observation confirmed that this simulated dynamics followed a FC model 
[c.f. ffTU]) ]. Finally, we observed that Rj(r) also exhibited a high degree of 
periodicity (Fig. Q] (e)), which was mostly due to the highly periodic activ- 
ity of Pxx{vi,vi,t). Accordingly, Sj(f) consisted of a sequence of harmonic 
frequencies separated by 1 Hz (Fig. 121(c)). 

4-2.2. Fully uncorrelated spatiotemporal dynamics 

In this dynamics, all the dipoles shared approximately the same nor- 
malized autocorrelation and total autocorrelation, p(v,v,r) ~ p{r) and 
Rj(v,v,r) ~ Rjij). Thus, this dynamics could be classified as ID [c.f. (TT41) 
and ffTB"]) ]. Furthermore, p{r) and Rj{t) were found to lack any periodicity 
and consisted of a single impulse located at r = 0, as illustrated in Fig. [[] 
(d) and (f). As opposed to the FC dynamics, P x (v) and P y (v) were approxi- 
mately the same. This observation is in agreement with the fragmented and 
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random nature of this simulated dynamics, in which dipoles were oriented 
with equal probability in every direction. 

The cross-correlation between dipoles decreased rapidly with the distance 
that separated them. Figured] (d) shows that p yy (vi,V2,r) consisted of a 
small but distinct spike, whereas p xx {v\, Vs,t) and p yy (vi,V3,r) were prac- 
tically zero. In other words, this dynamics exhibited highly localized spa- 
tiotemporal correlations. Accordingly, although this dynamics was approx- 
imately FU, it could be more accurately described as partially correlated 
(PC). We estimated the region whose dipoles showed a high degree of corre- 
lation with the dipole at v\, and found it corresponded approximately to A\. 
Interestingly, p xx (vi,V2,r) — 0, and the maximum oi p yy (vi,v 2 ,r) occurred 
at r = 0, as seen in Fig. [TJ (d). In words, only the y components of dipoles at 
v i and t>2 were correlated and they showed frequent simultaneous activities. 
The explanation for this observation is the following. Dipoles at v\ and v 2 
shared the same x coordinate and they were simultaneously activated when- 
ever local excitation wavefronts traveled along the y direction, which favored 
the generation of y dipole components. As a consequence, the y components 
of cardiac dipoles at vi and v 2 were often activated simultaneously and hence, 
p yy (vi, i>2, r) showed its maximum at r = 0. 

[Figure 2 about here.] 

[Figure 3 about here.] 

4-3. Analysis of the TDDF during FC dynamics 

Following ( 126]) . the spectrum of cardiac signals measured by pulse MSD 
during FC dynamics can be expressed as the product of the total spectrum by 
the Fourier Transform of the TDDF. We used our simulations to compare the 
spectrum of synthesized cardiac signals during FC dynamics, S Cn (f), and the 
spectrum predicted by (1261) . Firstly, we estimated spectra S Cn (f) by applying 
to c n (t) Welch's method [23j, with a 2048-samples Hamming window and 50% 
overlap. This configuration resulted in a high spectral resolution suitable to 
identifying the harmonic structure of c n (t). We subsequently built the time- 
delay function (a„(v,w) for each A n , by computing the difference between 
the excitation times of each pair of cells within A n . Then, we obtained 
the TDDF -Fa„(C) from CA n (v,w), and we calculated its Fourier Transform, 
SA n (f)- Finally, we obtained the spectrum predicted by (J2U|) by estimating 
the total spectrum of a single dipole Sj(f) and multiplying it by S'^ n (/). In 
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all cases, the spectrum obtained by virtue of ( 12 6 p matched the spectrum of 
synthesized cardiac signals. 

We illustrate these results with pulse MSD A 5 and A 8 . Figure[2](a) shows, 
on an isochronal map of cell excitations, the 21x21 pulse corresponding to 
A$, and the 81 x 81 pulse corresponding to Ag- The excitation times of 
four isochrones are also indicated, namely t = 0.01,0.05,0.07 and 0.11 s. 
According to this map, the earliest excitation in pulse A 5 (A 8 ) was t e = 0.05 
s {t e = 0.01 s), whereas the latest excitation was t\ = 0.07 s (ti = 0.11 s). 
Thus, the maximum time delay between excitations in pulse A§ (Ag) was 
(max = h -t e = 0.02 s (( max = t t - t e = 0.1 s). Also, since the time delay 
function satisfies ((v,w) = —((w,v), the minimum time delay in A 5 (A 8 ) 
was ( min = —0.02 s ((max = —0.1 s) and thus, the possible values of time 
delays within A 5 (A 8 ) lied in the range |C| < 0.02 s (|C| < 0.1 s). 

The generated TDDF, F A5 (() and ^(C), are shown in Fig. [2] (b). 
Because of the property ((v,w) = —((w,v), they were symmetric around 
zero. In addition to this, they exhibited a triangular shape and Fa 5 (() = 
(F As (() = 0) for time delays |C| > 0.02 s (|C| > 0.1 s). In all cases, we found 
that Fa„ (C) had a triangular shape, which is a consequence of the constant 
speed of the excitation wave. We also observed that the smaller A n , the nar- 
rower Fa„{(), since the maximum difference of time delays was also smaller. 
Accordingly, S^ n (f) had a square-sine shape and the smaller A n the broader 
^An(f)- As illustrated in Fig. [2] (c), Sj(f), S C5 (f) and S cg (f) shared the 
same harmonic structure, which consisted of harmonic frequencies located 
at multiples of 1 Hz. This is in agreement with the fundamental periodicity 
imposed by the simulated stimulation rate. Figure [2] (c) also shows that the 
envelopes of S Cg (f) and S C5 (f) conformed respectively to 5 , A g (/) and SA 5 {f)- 

4-4- Analysis of the spectral effects of lead systems 

According to (|26|) and (1281) . the spectral effects of pulse MSD are qualita- 
tively different for FC and ID-FU dynamics. In order to explore our analytical 
results, we compared the spectra obtained by pulse MSD A n during the sim- 
ulated FC and FU dynamics. We estimated the spectrum of each synthesized 
cardiac signal based on Welch's method, using a 512-samples Hamming win- 
dow and 50% overlap. This configuration resulted in a low spectral resolution 
suitable to study the spectral envelope. Based on the estimated spectra, we 
computed the 95% power bandwidth of synthesized signals during FC and 
FU dynamics (Tabled]). 
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During FC dynamics, we observed that as the size of A n increased, the 
bandwidth of synthesized cardiac signals decreased accordingly (Fig. El (a)). 
In order to assess quantitatively the effects of the LEV on the spectrum of 
cardiac signals, we compared the 95% power bandwidth of synthesized signals 
against the LEV (Tabled]). This comparison shows that during FC dynamics, 
the larger the LEV, the narrower the bandwidth. By contrast, during FU dy- 
namics the spectrum remained approximately unchanged across the different 
leads, except for the smallest pulses, Ai and A 2 (Fig. El (b)). Accordingly, 
the 95% power bandwidth remained approximately unchanged and discrep- 
ancies were observed only for the smallest LEV (Table [1]). The changes in 
the spectrum of cardiac signals synthesized by pulse MSD with small LEV 
during FU dynamics could be ascribed to the previously identified, localized 
spatiotemporal correlations. Within small regions of the cardiac source, it 
could no longer be assumed that the activities of the dipoles were uncor- 
rected and hence, signals measured by pulse MSD whose LEV was small 
enough, would exhibit a spectrum dependency on the LEV, characteristic of 
correlated sources. In our simulations, the estimated volume of correlation 
during FU dynamics was close to Ai and so, the effects of local correlations 
manifested more clearly in A\ and Ai. 

5. Conclusions and discussion 

In this paper we have presented a mathematical formalism for investigat- 
ing the spectrum of cardiac signals. Our formalism reveals that the spectrum 
of cardiac signals can be expressed in terms of the spectrum of the underlying 
cardiac source and the MSD of the lead system [c.f. ffT3l) ]. Two main con- 
clusions can be drawn from ffT3|) . Firstly, the information contained in the 
spectrum of cardiac signals is limited to the second-order characteristics of 
cardiac sources. And secondly, for the same underlying spatiotemporal dy- 
namics, different lead systems will, in general, produce cardiac signals with 
different spectra. Although the distortion introduced by the lead system fol- 
lows in a natural way from the linearity of the measurement (pQ), to the best 
of our knowledge this has not been formulated before. 

To gain understanding of the spectrum of cardiac signals, we have further 
analyzed (TT3|) for specific models of MSD and cardiac sources. Our analytical 
results show that the spectral effects of pulse lead systems is qualitatively 
different during FC and FU dynamics. During FC dynamics the spectrum 
of cardiac signals depends on the lead SR, which we quantified with the 
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LEV [c.f. (E)])]. Specifically, the spectral envelope of cardiac signals from 
local measurements (small LEV) is broad, whereas the spectral envelope 
from global measurements (large LEV) is narrow. By contrast, during FU 
dynamics the spectrum of cardiac signals is a volume average [c.f. f|27|) ] and 
hence, during ID-FU dynamics it remains the same irrespective of the LEV 
[c.f. (|28p ]. We have further explored the spectral manifestation of FC and 
FU dynamics in a simulation environment and our simulation results are in 
agreement with the analytical ones. 

Based on our simulation experiments, we could also explore the spectral 
manifestation of PC sources. In the analysis of the spatiotemporal char- 
acteristics of simulated FU dynamics, we observed localized spatiotemporal 
correlations which allowed us to more accurately describe the simulated dy- 
namics as PC. Also, when analyzing the spectrum of cardiac signals synthe- 
sized from this dynamics, we observed spectral distortions for small LEV. 
These observations could be explained as follows. For large LEV, lead sys- 
tems measure uncorrelated dipoles, and no spectral distortions are observed. 
However, for small LEV lead systems focus within regions where the dipoles 
are in fact correlated and accordingly, some distortions on the spectrum of 
cardiac signals are expected. Our observations encourage us to speculate that 
in general, irrespective of the type of MSD, the lead SR is a major factor in 
determining the spectrum of cardiac signals. We have derived analytically 
this result for idealized pulse MSD during FC and FU dynamics, but not for 
generic MSD nor for PC dynamics. We can speculate that in those scenarios 
where the LEV is larger than the volume of correlation, PC dynamics can be 
treated as FU, whereas whenever the LEV is within the volume of correla- 
tion, PC dynamics can be treated as FC. Nevertheless it remains to further 
study the spectrum of cardiac signals for generic lead systems and generic 
spatiotemporal dynamics. The formalism that we have presented constitutes 
a convenient mathematical framework for investigating the spectrum of car- 
diac signals measured during generic spatiotemporal dynamics. Within our 
formalism, a meaningful definition of the degree of spatiotemporal correla- 
tion and the volume of correlation could be proposed, and by using (TT5j) its 
effects on the spectrum of cardiac signals could be thoroughly investigated. 

The study that we have presented has practical implications on the analy- 
sis of the spectrum of cardiac signals. We have shown that during correlated 
dynamics leads can distort the spectral envelope, whereas the harmonic sepa- 
ration remains unchanged. Therefore, spectral features which depend on the 
spectral envelope such as the peak frequency (g,lZ|, the mean frequency 24J, 
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and the organization index |4| , can be affected by the characteristics of the 
lead system, whereas those features that only depend on the harmonic sepa- 
ration, notably the DF [8|, |9j, [l0|, llll, ll2[ , are insensitive to the characteristics 
of the lead system. This conclusion supports, from a theoretical standpoint, 
recent experimental studies which assessed the robustness of the DF. In [25[ 
it was reported that the values of the DF obtained from unipolar and bipo- 
lar cardiac signals correlate highly with the DF values obtained from optical 
signals. Other studies have compared the DF values obtained from cardiac 
signals measured by contact and non-contact intracardiac mapping systems, 
and have found good ag reement between the DF values obtained from each 



mapping system [17|, |18|. Finally, a high correlation between DF values from 
intracardiac EGM and surface ECG signals was found in 26(. In summary, 
these experimental studies indicate in agreement with our theoretical results, 
that the DF is a robust spectral feature of cardiac signals. Therefore, there 
is both experimental and theoretical evidence of the robustness of the DF, 
which pave the ground for developing novel DF mapping systems for di- 
agnosing and managing of fibrillation, such as non-contact, intracardiac DF 
mapping systems 17|,ll8| and non-invasive, surface DF mapping systems 19 . 



Appendix A. Derivation of the Correlation and the Spectrum of 
Fully Correlated Sources 

By setting v = w in ([3]) and substituting (ITS]) into it, it can be proved that 
the autocorrelations of any two dipoles 3(v,t) and J(w,t) in a FC source, 
are always identical. 

p(w,w,r) = (J' T (w,t + T)J'(w,t)) t 



(J' 1 (v, t + r - C(v, w))3'(v, t - C(v, w))) t 



r'T7„. V 



(J'>,i' + T)J'(i;,0>< 

p(v,v,r) 



(A.l) 



where t' = t — ((v,w). Therefore, FC sources are by construction also ID, 

p(v,v,r) = p{r) and er(v,v,f) = cr(f). By substituting ffT8|) into ([3]), 
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p(v, w, t) can be expressed in terms of pij): 

P (v jW ,t) = (J' T (v,t + r)3'(w,t)} t 

= (3 /T (v,t + r)J'(v,t-av,w))) t 

= (J' T (v,t + T)J'(v,t)) t *6(T-av,w)) 

= p(r) * S(r — C(v,w)) 

= p(T-((v,w)). (A.2) 

In words, p(v,w,r) is a ((v, u>)-delayed version of pij). As a consequence, 
by applying the Fourier Transform to (1A.2J) . the spectrum of a FC source can 
be expressed as: 

<r(v,w t f) = F\p(v,w,t)] 

= F[p(t)*6(t-((v,w))} 

= F[p(t)]T[6(t-((v,w))} 

= <r{f)exp[-j2nfC{v,w)]. (A.3) 

Appendix B. Connecting the Spectrum of Cardiac Signals to the 
Spectrum of Cardiac Sources 

By substituting (TjQ) into (jHJ), we can express c'(t) in terms of J'^, t), 

c'(t) = I L T (v) J(v, t)dv - ( f L T (v) J(v, t)dv) t 
Jv Jv 

= J L T (v)J(v,t)dv - I L T (v)J T (v)dv 
Jv Jv 

= f L T (v)(J(v,t)-3(v))dv 
Jv 

L T (v)3'(v,t)dv. (B.l) 



v 



Then, by substituting (IB. 1|) into (fTU|) . we can write the autocorrelation of 



a cardiac signal R c {t) in terms of the autocorrelation of the cardiac source, 
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p(v,w,r): 

/?,(?) (' / L T (v)3'(v,t + r)dv J L T ( y w)J'(w,t)dw) t 

Jv 



l c V 



/ L T (v) (3'{v, t + t)3' t (w, t)) t L(w)dvdw 

JVxV 

L T (v)p{v,w,T)L(w)dvdw. (B.2) 



VxV 



Finally, following (TTTj) . we apply the Fourier Transform to (IB. 2p and derive 
the relationship between the power spectrum of a cardiac signal S c (f) and 
the spectrum of the cardiac source, (t(v, w, /): 

S c (f) = J~ \ \ Li T (v)p(v,w,r)'L(w)dvdw> 

{JvxV ) 

= I \j T {v)J r {p{v,w,r)}\j{w)dvdw 

JVxV 

L T (v)a(v,wJ)L(w)dvdw. (B.3) 



VxV 
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List of Figures 

1 Analysis of simulated fully correlated (FC) and fully uncor- 
rected (FU) dynamics. (See text for explanation.) Voltage 
maps during (a) FC and (b) FU dynamics, showing source 
locations v%, v 2 and v 3 . Normalized correlations P xx (vi,w,t) 
and p yy (vi,w,r), where w = Vx,V2,Vs, during (c) FC and (d) 
FU dynamics. Normalized total autocorrelation Rj(t) during 



(e) FC and (f) FU dynamics |24 

The envelope of cardiac signals measured during fully corre- 
lated (FC) dynamics conforms to the Fourier Transform of the 
time-delay density function (TDDF). (a) An isochronal map 
of cell excitations during simulated FC dynamics and the re- 
gions covered by pulses A 5 and A$ are shown, (b) Based on 
the time delay between cell excitations, the TDDF correspond- 
ing to pulses from leads A 5 and A 8 , Fa 5 (() and Fa 8 ((), are 
generated, (c) As predicted by (|26|) . the multiplication of the 
total spectrum Sj(f) by the Fourier Transforms of the TDDF 
SA s (f) and SA 8 {f), matches the spectra of signals generated 



by pulses A 5 and A 8 , S CB (f) and S cs (f) |25 

The spectral effects of the characteristics of pulse MSD are 
qualitatively different during FC and FU dynamics, (a) Dur- 
ing FC dynamics, the smaller the pulse, the broader the band- 
width of synthesized leads, (b) During FU dynamics, spectra 
remain approximately unchanged across different leads except 
the smallest ones, A\ and A 2 . This behavior can be explained 
by the existence of localized spatiotemporal correlations in the 



simulated FU dynamics |26 
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Figure 1: Analysis of simulated fully correlated (FC) and fully uncorrelated (FU) dy- 
namics. (See text for explanation.) Voltage maps during (a) FC and (b) FU dynam- 
ics, showing source locations v\, v^ and V3. Normalized correlations /} xx (vi,w,t) and 
p yy (vi,w,T), where w = Vi, 1)2,1)3, during (c) FC and (d) FU dynamics. Normalized total 
autocorrelation Rj{t) during (e) FC and (f) FU dynamics. 
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Figure 2: The envelope of cardiac signals measured during fully correlated (FC) dynamics 
conforms to the Fourier Transform of the time-delay density function (TDDF). (a) An 
isochronal map of cell excitations during simulated FC dynamics and the regions covered 
by pulses A$ and As are shown, (b) Based on the time delay between cell excitations, the 
TDDF corresponding to pulses from leads A$ and As, Fa 5 (C) and Fa 8 (£), are generated. 
(c) As predicted by (|2l)l) . the multiplication of the total spectrum Sj(f) by the Fourier 
Transforms of the TDDF Sa 5 {J) and SU S (/), matches the spectra of signals generated by 
pulses A 5 and A s , S C5 (f) and S Cg (f). 
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Figure 3: The spectral effects of the characteristics of pulse MSD are qualitatively dif- 
ferent during FC and FU dynamics, (a) During FC dynamics, the smaller the pulse, the 
broader the bandwidth of synthesized leads, (b) During FU dynamics, spectra remain 
approximately unchanged across different leads except the smallest ones, A\ and A^. This 
behavior can be explained by the existence of localized spatiotcmporal correlations in the 
simulated FU dynamics. 
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Tabic 1: Definition and lead equivalent volume (LEV) of the measurement sensitivity 
distributions used in the simulations and bandwidth of the corresponding simulated signals 
during fully correlated (FC) and fully uncorrclated (FU) dynamics 



Lead 


Dimensions 


LEV (%) 


BW FC (Hz) 


BW FU (Hz) 


A 


5x5 


0.2 


114 


79 


A 2 


9x9 


0.8 


77 


63 


A 3 


13 x 13 


1.6 


85 


61 


A, 


17 x 17 


2.8 


70 


60 


A b 


21 x 21 


4.3 


60 


58 


A e 


41 x 41 


16.4 


35 


54 


A? 


61 x 61 


36.4 


28 


54 


A 8 


81 x 81 


64.3 


24 
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